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' Exclusive double charmonium meson production in e+e" annihilation at energy y/s = 10.6 GeV remains very 
^ ] interesting problem for theoretical investigations. NRQCD that is often used to predict the cross section of such 
processes fails to archive agreement with Belle and BaBar experiments. The predictions for the cross sections obtained 
in the framework of NRQCD0 are by an order of magnitude less than the experimental results 13, 

An interesting approach to the problem was proposed in papers^ 6J. In these papers the process e^e~ iprjc 
at energy -^/s — 10.6 GeV was considered in the framework of light cone formalism. In paper p( it was shown that 
taking different light cone wave functions one can considerably enhance NRQCD prediction for the cross section of 
the process e+e^ tpTjc- Using physical model for the light cone wave functions of ip, rjc mesons the authors of paper 
, received the prediction for the cross section of this process that agrees with experimental result obtained at Belle. 

• In addition to the agreement between the theoretical prediction and the experimental results light cone approach 
Cn| allows one to understand that charmonium mesons wave functions are two wide to be considered in the framework of 

NRQCD in such processes. 

Further progress in the understanding of exclusive double charmonium mesons production is connected 
. with paperQ. In this paper it was shown that contrary to NRQCD the study of the processes e+e^ 
' ^{lS)r]c{2S),'if{2S)ric{lS)j_^{2S)ric{2S) in the framework of light cone formalism leads to a good agreement with 
Ph. Belle and Babar results 0,111 • 

• So light cone formalism better than NRQCD predicts the cross sections e+e^ 

(]J ' J/*?7c(l'5'), J/*?7c(2S'), *(25')?7c(lS'), *(25)?7c(2S'). In addition to these reactions Belle and BaBar have mea- 
sured the cross section of the processes e~^e^ J/'i>XcO,'^{2S)xcQ- The agreement of NRQCD prediction and 
^ ■ experimental results is also poor for this reactions. The aim of this paper is the application of light cone formalism 
to these reactions. 

This paper is organized as follows. In the next section light cone wave function of J/'I', ^'(25') and XcO mesons 
will be considered. In the third section the expression for the amplitudes of the processes under consideration will be 
derived. Numerical analysis of the cross sections e~^e~ — > J/^XcOj ^(2>5')xcO will be presented in the forth section. 
Finally the results of the calculation will be discussed in last section. 



II. LIGHT CONE WAVE FUNCTIONS OF J/*, *(25') AND XcO MESONS. 

To calculate the amplitude of the process e+e^ — > l^Xco[^ = J/'^,'^{2S)] one needs to know light cone wave 
functions of the final charmonium mesons. The functions for J/^E* and ^'(25') mesons are defined as follows^ 
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The dependence of the hght cone wave functions on the scale fi is very slow and it will not be considered in the full 
form. Only renormalization factors of the corresponding local currents will be regarded. The constants /* , can be 
determined from QCD equations of motion 
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where Mq = Mq^{^i — Mq^). The factors Zp.Zt, Z„i can be written in the form 
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where Cp = 4/3, bo = 25/3. 

For the light cone wave functions of J/^* meson the model proposed in 1^ will be used: 



X1X2 



[1- 4X1X2(1 

For the light cone wave functions of ^'(25') meson the model proposed in 7J will be used: 



(1 - V^)xiX2 



X1X2 



[1 - 4x1x2(1 - v^)] / [1 - 4x1x2(1 - V^)] 
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Here w is a characteristic speed of quark-antiquark pair in meson, ci is the coefficient fixed by the wave function 
normalization condition J dx0i(x,f^) — 1. The constant /3 is fixed by the condition that zero of the wave function of 
25 state meson 7] must coincide with zero obtained from the solution of Schrodinger equation with BuchmuUer-Tye 
potentially ■ 0° are the asymptotic expressions of the wave functions 



Vk(x) = Vl(x) = Vrix) = (pasyix) = 6X1X2, 



(6) 



y^(x) = -[l + (xi-X2)2]. 



(7) 



It should be noted here that the expressions for light cone wave functions Q,© link different limits: quark-antiquark 
pair in meson being at rest v ^ and very light quark v ^ 1. In the former limit one obviously gets ~ S{x — 1/2), 
the later one leads to the asymptotic function ~ (j>°'. 

Let us consider the light cone wave functions of Xco meson: 



(Xco(p)|Q(^)g(-^)|0> ^fsj dye'p^^y^-y^^ssiy), 

(Xco(p)|Q(^)7mQ(-^)|0) = dye^P^^y^-y^'^Sv{y) + fl''^z^J dye^P^^y^-y^^S^^\y), 

{Xco{p)\Q{z)a^,Q{-z)\0)^fT{p^.z,-p,z^) [ dye'P<y^-y-^ Sriy). 



(8) 

The functions (pi — Ss{y), Sy \y), Sxiy) are normalized as follow J dy<l)i = 1, the normalization condition for Sviy) 
is / dy{yi — y2)Sv(y) — 1- Using QCD equation of motion the constants fs-, /t, fy^ can be related to fy^: 
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(9) 
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The constant fs and consequently /y ^ can be expressed through matrix element (Xco(j')|Q(0)Q(0)|0) found in 
As the result we get 



where i?'(0) is derivative of the wave function of heavy quarkonium at the origin, mq is the mass of c quark in 
potential model, < Oi >p is well known from NRQCD matrix element 0|. It should be noted that formulae ® 
were found in paper ^| but the expression for the constant fy ^ differs from (|10|l and agrees with [T^ . 

Although there are four light cone Xco meson wave functions only Ss{y) and Syiy) contribute to the amplitude 
of the processes at accuracy considered in our paper, the others give power correction to the result. Below we will 
consider only Ss{y) and Sviu)- 

As in the case of vector mesons we are not going to regard full /i-dependence of light cone wave functions. Only 
overall renormalization factor of corresponding operator will be considered. The rcnormalization factor for light 
cone wave function Ss{y) equals to Zp(fj,). To find renormalization factor for the function Sv{y) it is worth noting 
that the corresponding operator (Xco(p) 10(0)7^(5(0) |0) equals zero. So the renormalization factor for the func- 
tion Sv{y) equals the renormalization factor of the first nonvanishing multiplicatively renormalized local operator 
{Xcoip)\Ql^.Cl^\zD/zd)Q\0). It equals to[ll| 



L«.(^q)J 



8Cp 
9bo 
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In order to get the expressions for the light cone wave functions Sv{y) and Ss{y) we will use the same procedure 
that was used in paper |^ to get expressions Q). 

First let us consider the leading twist wave function of Xco meson Sviy)- There are three wave functions of 2P 
meson with different orbital momentum projection Lz in the direction of meson motion(m = ±1, 0). But only m — 
projection gives contribution to the leading twist light cone wave function 11] . This wave function can be approximated 
by Coulomb wave function of 2P state: 

where v is the characteristic relative velocity. Substituting^^ 

k^^k±, fc.- (2/1 -2/2)^, Mo2 = I^i±^ (13) 

2 yiy2 

and than carrying out the integration 

' d'k^^,iy,k^) (14) 
one gets the expression for the light cone wave function: 



(P- Sy i yiy2-7z — -j ^ 

[ (1 - 4yiy2(l - -y^))^ 



(15) 



In last formula the expression Sy^ = 1/12/2(2/1 ~ 2/2) is asymptotic form of the light cone wave function Sviy)- To get 
final expression for the wave function we will modify H15() similar to wave functions of 15 and 25' mesons (@J, (0: 

Sviy) = cv{v^)St^^{y)!y,y2 ~ ^^'^^ ' ""'L ] (16) 

where cv(w^) is the coefficient which is fixed by the wave function normalization J dySviy)iyi — 2/2) — 1- For the 
Ss iy) wave function the following expression will be taken 

l~2ymil-v^) 



Ssiy) = cs(«^)5r(2/) 2/12/2 ,^ ."^7. , (17) 

I (1 - 42/12/2(1 -w^))"' J 



where asymptotic form of the wave function 5^*(j/) — 1, cgiv'^) is the coefficient which is fixed by the wave function 
normalization J dySsiy) — 1- 
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III. THE CALCULATION OF THE AMPLITUDE OF THE PROCESSES e+e" ^ V"(pi, ei)xco(p2). 
Leading asymptotic behavior of the matrix element ^V{pi,ei), XcaiP2)\J^''\0) can be derived from formula [l^ 

(A.f(pi,Ai)Af(p2,A2)|j;'|0) ^ (^i=y ' , (18) 

where J^*^' is the electromagnetic current. For the processes under consideration we have M{pi,Xi) = 
V{pi, Xi), M{p2, X2) — Xco{P2)- Obviously the helicity A2 equals zero. As to the vector meson V{pi,Xi) the leading 
contribution is given by the helicity Ai = 0. So the asymptotic behavior of the amplitude is 

(F(Pi,Ai),^(p2)K|0) ^ ^ (19) 

and the asymptotic behavior of the cross section (T(e+e" — > VS) is ~ 1/s'^. Unfortunately one can show that it is 
not possible to disregard NLO contribution in 1/s expansion. To see this let us consider NRQCD result for the cross 
section of the process e+e^ J/'(p{lS)xcO obtained in paper 0. Cross section of this process can be written as 

"""^ 2 2 2 17 2 f, 'ofvf's /on^ 

a = -T-a a.q^For a/I-H-^, (20) 

where = 16ml/ s and Fq = 2{l8r'^ - Ir'^f +r'^{A+ lOr^ - Sr^^. Let us substitute s 10.6^^ and expand the 
above formula in 1/^ series (numerical inputs are the same as in Q). We get 

0.23 2.91 0.89 

- = ^ + — - — + (21) 

Thus one sees that in the framework of NRQCD NLO correction at energy ^/s = 10.6 GeV is by an order of magnitude 
larger than the leading one. So one can suppose that NLO term in 1/s expansion gives considerable contribution and 
must be regarded in our analysis. The same is true for the process e"*"e~ ip(2S)xcQ- 

Two diagrams that contribute to the processes e~^e~ — *■ V{pi, Xi)xcoiP2) are shown in Fig.l. and the other two 
can be obtained from them by charge conjugation. Having the expressions for the light cone wave functions |(SJ), 
((T7|) it is not difhcult to obtain the matrix element {V{pi, Ai), Xco(P2)|«/^|0): 

{V{pi,Xi)Xco{P2)\J^\0) - .9i(p5' -p^)(eA,P2) + .92((eAiP2)K - <(piP2)) (22) 
where the formfactors gi and 52 are: 



- 2^^M.M, I' I' ,y,aM (-16^.(^)™^^ ) (23) 



MvM^ / dxi / dyias(k)[2Zy{k) — (I - Z,n[k)Zt[k)——^)— 



— 2 

l6ZmZ.{k)ZUk)-^-^^^^^ - 8Z.{k)^^-^^iVUx) 2V^{x) - + 



32^.(fc)^*(^)^^e^^/), (24) 



s{x)d{x,y) Ml 

where d{x, y), s{x), sly) are dimensionless quark and gluon propagators defined as follows 

fc2 / s\ _ [Z^ik)MQ' 



d{x,y) = ^= xi + — h;i + — , ^ , (25) 



s 



V yiy2s J \ X1X2S J 

where fc is a momentum of virtual gluon, cr is a characteristic momentum of virtual quark in the diagrams shown in 
Fig.l, the constant / = -3M^^J2'Mq. 
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It is interesting to note that there are two contributions to the formfactors gi and 52 factored in formulae H23|l . 
(|24|l . The first is NRQCD result for the formfactors that does not regard internal motion of quark-antiquark pair 

inside mesons and it is proportional to ~ fy ^ fv ^ R{0)R'{0). The second contribution regards internal motion and 
it is proportional to the integrals / dxidyi in expressions (|2SJ, iP^ - 

In the limit v ^ the mesons V and Xco have equal masses M — My = Af^ and formulae (|23|l . (|24ll must reproduce 
NRQCD result for the process e+e^ J/'^XcO- 

Really if one takes the limit v ^ m expressions H23(l . 124(1 one gets 



s 



So one sees that at accuracy considered in our paper the expressions ()27|l and 1(28(1 coincide. 
The cross section of the processes e~^e~ Vxco is given by the formula: 

<»' 

where qc is the charge of c-quark, pi is the momentum of vector meson in final mesons' center mass frame, for the 
electromagnetic current of the form ((22(1 the function F is given by the formula 



17 9lQo I S2I0 1 6 I -'^ 2 11 r2 4 4 2A/r2 4 I 1 27i,r2 4 ^ 71 4 

IglM^Mlql + ^g^g^M^M^ql + ^ff^M^M^ + ^^'^yA^v + lgig2M^M^ 



2„4 



5^54 2g2s3 2^3 2glMls^ 1 , , 

^1 ^1 ' ^ .913253 + 0(l/s2), (30) 



where ql — s — My — M. 



6M2 3 3 3Af2 3 



IV. NUMERICAL CALCULATION. 



In the numerical analysis the following parameters will be used: 



Mc = 1.2GeV, 
^(15), /y = 0.41GeV, 

tp{2S), fv = 0.28GeV. (31) 



The values fv were obtained from decay width 



As to the value of the constant fy'' it is related to the value of the constant fs = (Xco(p)|Q(0)Q(0)|0) found in 
using QCD sum rules: 

/^i) = 0.084.GeV (33) 
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H\H2 


CySaBar X Bt H2-tcharged>2{^)^ 


O Bella X -BrH2^charsed>2(fb)f3] 


CTLO (fb) 


O-iVflQCD (fb)[2j 


^(IS')xco 


10.3±2.5l:i;8 


6.4 ± 1.7 ± 1.0 


14.4 ifa'a 


2.3 


^(2S)xco 




12.5 ±3.8 ±3.1 


7.8 f:i 


1.0 



TABLE I: The second and third column contain experimental result. The results ol our calculation is presented in the forth 
column. The last column contains NRQCD results. 



For as {jj) one loop result will be used 



47r 

= bolog(MVA^) ' ^^'^ 



with A = 200MeV. The other parameters needed for the numerical analysis is the width of wave function . It will 
be taken from potential model jg: 

J/^ = 0.23, 
Xco i^' = 0.25, 

%p{2S) = 0.29. (35) 

Now let us consider the formula of the cross section (|29|l . It is seen from (|30|l that leading order contribution(LO) to 
the cross section is given by the first term ~ gl and the rest of the formula H30I) is the power correction to the LO. 
In our paper the power correction to the formfactor gi is beyond the accuracy of our calculation. But it is well seen 
from H30|l that in order to get the cross section up to 0{l/s^) one must know 1/s correction to the formfactor gi. 
Fortunately LO term ~ g\ gives negligible contribution to the cross section and varying g\ in a reasonable region the 
cross section is changed by few percent. Sure we do not pretend to the accuracy about few percent in our calculation. 

In our paper we use the model for the light cone wave function defined by equations Q), (|SJ|, (|16l) and H17|) 
with the widths (|35|1 . To estimate the size of this uncertainty results from this model in addition to the widths 
(|35|l the calculation of the cross sections with shifted widths(about 10%) is done. The results of the calculation 
is presented in Table 1. The second and the third columns contain experimental results measured at BaBar and 
Belle experiments. In the fourth column the results of this paper are presented. Central values of the cross sections 
correspond to the light cone wave functions with widths (|35|l . The upper values of the cross section correspond 
to the widths 0.26( J/^*), 0.32(\E'(2S')), 0.28(xco)- The lower values of the cross section correspond to the widths 
0.2( J/^f), 0.26(5'(25)), 0.22(xco)- In order to compare the result with NRQCD predictions for the processes under 
consideration the fifth column contains the predictions in the framework of this model. 

From Table 1 one sees that the predictions of the cross section of the processes e+e" J/^XcO, *(2S')xcO in 
the framework of light cone are larger than NRQCD predictions and the agreement with the experimental results is 
better. As it was noted in |2j the difference of the NRQCD and light cone prediction for the cross sections can be 
attributed to the fact that at leading approximation NRQCD does not regard the motion inside final mesons. So 
NRQCD predictions for the cross section of the processes under consideration are unreliable. 

In addition to the uncertainties described above and uncertainty due to the unknown size of radiative QCD correction 
one can suppose that there is very important 1/s correction. The size of this correction can be estimated from formula 
(|21|l for NRQCD result of the cross section of the process e+e" J/'^XcO- It is seen from H21f) that in the framework 
of NRQCD 0{l/s^) contribution changes the value of the cross section by 30 %. Moreover 0{l/s^) correction diminish 
the value of the cross section. It was noted above that to get light cone result one must multiply NRQCD by a factor 
regarding internal mesons' motion. If one suppose further that these factors for 0(l/s^) and 0(l/s^) contribution 
are of the same order than one can claim that in the framework of light cone formalism 0(l/s^) contribution diminish 
the size of the cross section by 30%. After including this correction the value of the cross sections can be estimated 
as cr(e+e- J/*Xco) ~ 9fb and cr(e+e- ^ ^'(2S')xco) ~ 5fb. 

It should be noted here that light cone prediction for the cross section of the process e+e^ '^{2S)xcO is almost 
twice less than Belle result. One can attribute the difference to any source of uncertainty described above. But 
another source the disagreement can arise from the higher fock state of XcO meson. Really it is known from NRQCD 
that color octet contribution of XcO meson is of the same order in relative velocity expansion as color singlet. Moreover 
it is known that color octet state gives NLO contribution to the amplitude, i.e. it changes main contribution to the 
cross section. So it would be interesting to estimate the size of this correction. 
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HiH2 


D /-Cu. \ r /1 1 

OBaBar X -DrHa ^charsed>2 (.IDJ 1^ 


O'Belle X -Brfl-2^char9ed>2(fb)[^ 


crLo(fb)[7] 


CTNRQCD (fb)[2j 


^(15)77,(15) 


17.6 ±2.812;? 


25.6 ±2.8 ±3.4 


26.7 


2.31 


^(2S)j7c(lS) 




16.3 ±4.6 ±3.9 


16.3 


0.96 


V'(lS)»7c(25) 


16.4 ±3.7^3:0 


16.5 ±3.0 ±2.4 


26.6 


0.96 


^(2S)»7c(2S) 




16.0 ±5.1 ±3.8 


14.5 


0.40 



TABLE II: The second and third column contain experimental result. The results of paper Q are presented in the forth column. 
The last column contains NRQCD results. 



DISCUSSION. 



In this paper the calculation of the cross sections of the processes e+e — > J/'^XcO and e+e 



*(25)XcO at 



energy y/s = 10.6GeV in the framework of light cone formalism has been carried out. It is shown that regarding the 
internal motion of mesons in the hard part of the amplitude in the framework of light cone results to the considerable 
enhancement of the cross sections in comparison to the NRQCD where internal motion is disregarded. So NRQCD is 
unreliable for the calculation of these cross sections. 

In addition to the processes e+e^ J/^'XcO, ^(2S')xcO Belle and BaBar experiments have measured the cross 
sections e+e" — > J/'^rjc, ^j2S)r]c, J /'^rjc{2S), {2S)rjc{2S) . In the framework of light cone formalism these reactions 
were considered in paper The results obtained in this paper are presented in Table 2. Comparing the results for 
the cross sections measured in Belle and BaBar with light cone and NRQCD predictions one can claim that despite 
a number of uncertainties the results obtained in the framework of light cone are in better agreement with Belle and 
Babar results than NRQCD predictions. 
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